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Abstract — Filter Bank-based Multicarrier Modulation (FBMC) using Offset Quadrature 
Amplitude Modulation (OQAM), known as FBMC/OQAM (or OFDM/OQAM), provides an 
attractive alternative to the conventional Cyclic Prefix-Based Orthogonal Frequency 
Division Multiplexing (CP-OFDM), especially in terms of increased robustness to frequency 
offset and Doppler spread, and high bandwidth efficiency. It suffers, however, from an 
inherent (intrinsic) imaginary Intercarrier/Inter-symbol interference that complicates signal 
processing tasks such as Channel Fstimation (CE). Recently, the so-called Interference 
Approximation Method (IAM) was proposed for preamble-based CE. It relies on the 
knowledge of the pilot's neighbourhood to approximate this interference and constructively 
exploit it in simplifying CE and improving its performance. The IAM preamble with nulls at 
the neigh boring time instants and extended version of it, which can provide significant 
improvement through an appropriate exploitation of the interfering symbols from 
neighbouring time instants that results in CE performance was recently reported. In this 
paper, we present IAM preamble design and apply it to identify the optimal IAM preamble 
sequence which results in a higher gain. Numerical results have verified the effectiveness of 
the theoretical framework and a gain of 1.24 dB against E-IAM-C. 

Index Terms— FBMC/OQAM, CE (Channel Estimation), IAM (Interference 
Approximation Method), E-IAM-C (Extended IAM-C), Pilots. 

L Introduction 

Orthogonal frequency division multiplexing (OFDM) has become quite well-liked in each wired and wireless 
communications [4, 2], in the main owing to its immunity to multipath attenuation, which allows a significant 
increase in the transmission rate [13]. using the cyclic prefix (CP) as a guard interval, OFDM manages to 
show a frequency selective channel into a set of parallel flat channels with independent noises. This greatly 
simplifies each the estimation of the channel and the recovery of the transmitted data at the receiver. 
However, these advantages come at the cost of an enhanced sensitivity to frequency offset 
and Doppler spread. This can be as a result of the fact that, although the subcarrier functions are perfectly 
localized in time, they suffer from spectral leak in the frequency domain and hence inter-carrier interference 
(ICI) results. Moreover, the inclusion of the CP entails a waste in transmitted power further as in spectral 
efficiency, which, in practical systems, can go up to twenty five percent [2]. 

Filter bank-based multicarrier modulation (FBMC) using offset quadrature AM (OQAM), known as 
FBMC/OQAM or OFDM/OQAM [8], provides an alternative to CP-OFDM that can mitigate these 
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drawbacks. FBMC/OQAM employs pulse shaping via an IFFT/FFT based efficient filter bank, and staggered 
OQAM symbols, i.e., real symbols at doubly the symbol rate of OFDM/QAM, are loaded on the subcarriers 
[12]. This allows for the pulses to be well localized in each the time and the frequency domains while still 
keeping most spectral efficiency. As a consequence, the system's robustness to frequency offsets 
and Doppler effects is enhanced and at the same time an enhanced spectral containment, for information 
measure sensitive applications (e.g., cognitive radio [1]), is offered. Moreover, 
FBMC/OQAM doesn't need the inclusion of a CP, which may result in even higher transmission rates [12]. 
However, the subcarrier functions are now only orthogonal in the real field, which 
means that there's perpetually an intrinsic imaginary interference among (neighboring) subcarriers and 
symbols [7]. As a consequence, the simplicity of the channel estimation (CE) in CP-OFDM is lost in 
FBMC/OQAM systems, because the real-valued pilots also are polluted" by this imaginary- valued 
contribution of their neighbours. A number of training schemes and associated channel estimation methods 
that take this under consideration have recently appeared in the literature, together with each preamble-based 
(e.g., [9,11]) and scattered pilots-based (e.g., [7,10,3]) ones. These can be roughly characterised as aiming at 
cancelling the unwanted interference or constructively exploiting it to boost the estimation performance. The 
latter approach depends on the fact that the interference to a pilot symbol is principally contributed to by its 
first-order neighbors, and hence it may be approximated if these neighbours also carry known symbols, a not 
uncommon situation within a preamble. This has been known as the Interference 
Approximation method (IAM). The idea is to combine the real-valued pilot with the interference 
approximation into a complex-valued pseudo-pilot and performing channel estimation as in CP- 
OFDM, beneath the idea of a locally constant channel frequency response. As the pseudo pilot divides the 
noise components, its magnitude should be as large as possible. This way, the interference turns into a benefit 
for the channel estimation task and hence the neighboring pilots should be therefore chosen on maximize 
its effect. 

A number of variants of IAM have been suggested, corresponding to different alternatives for the preamble 
structure [18]. Those that have obtained most of the attention comprise three FBMC/OQAM symbols (i.e., 

I. 5 complex OFDM/QAM symbols), with the edge ones being all zeros. This is often to safeguard the 
center vector of pilots from being interfered by the unknown knowledge components of the previous and 
current frames the center pilot symbols are then chosen therefore on maximize the interference contributions 
from neighboring subcarriers. Recently, such a preamble structure with imaginary pilots was suggested in 
[16] and individually in [6] and shown to be the optimum alternative inside this kind of preambles. It was 
called IAM-C to signify the presence of complex (imaginary) pilots, in analogy with IAM-R [9], the 
optimum real-valued preamble. IAM-C advances upon IAM-I [11], a preceding IAM variant furthermore 
employing invented pilots, by maximizing the magnitude of all pseudo-pilots rather than of only one third of 
them as in IAM-I. 

The possible profits from a preamble that has nonzero pilots at its side FBMC/OQAM symbols as well. 
Recently, a preamble structure [17] of this kind is proposed that exploits the symmetry of the interference 
weights and their relation values to get pseudo-pilots of a magnitude larger than in IAM-C. Recently, the 
expands IAM-C [17] to incorporate edge symbols as well is proposed, it will be called expanded IAM-C (E- 
IAM-C). 

Motivated by the standard of IAM-C method and awareness of the nature of the three IAM methods in [17, 
11], in this paper we have formulated a general theoretical structure for IAM preamble design. As a outcome 
the optimal preamble IAM sequence has been recognised which results in a higher performance gain. The 
effectiveness of the theoretical framework and the superiority of proposed preamble has thereafter been 
verified by numerical simulation. 

The rest of the paper is organized as pursues: In part II, the discrete -time baseband equivalent model for the 
FBMC/OQAM scheme is recounted. Part III reconsiders the three major preamble structures that encompass 
edge none guard symbols, IAM-R, IAM-C and E-IAM-C. The suggested optimal preamble is developed and 
considered in part IV. Replication outcomes are described in part V. Part VI concludes the paper. 

II. System Model 

The (QAM or OQAM based) OFDM modulator output is transmitted through a channel of length L h , which 
is assumed to be invariant in the duration of the preamble. At the receiver front-end, noise is added, which is 
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assumed Gaussian with zero mean and variance cr 2 . For CP-OFDM, the classical Configuration is assumed 
[13]. The discrete -time signal at the output of an FBMC/OQAM synthesis filter bank (SFB) is given by [12] 

M-l 
m = n 

Where d mn are real OQAM symbols, and ^ _ ^ ^ _ n g y lT*( ! 2~ ) g ;»>„,.„ with g being the 
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real symmetric prototype filter impulse response (assumed here of unit energy) of length L g , M being the 
even number of subcarriers, and ^ = (p fj +—{m + n)moArc whereto can be arbitrarily chosen [12]. In 

this paper, and without loss of generality, (p m n is defined as ^ m + n ^ — - mn7l as in [12]. The filter g is 
usually designed to have length jr = KM ,with K being the overlapping factor. Here the transmitted 
symbols d m n are real-valued with symbol duration T and inter carrier spacing V respectively. Two kinds 
of realisations of OFDM/OQAM are of practical interest. One can either set V = F , T = T 1 2 as in [14], 
[15] or set T Q = T , v = F I 2 as in [14]. We use the former approach. The double subscript (,) m n denotes 
the (m, ri) * frequency time (FT) point. Thus, m is the subcarrier index and n the OQAM symbol time index. 
The pulse g is designed so that the associated subcarrier functions g m n are orthogonal in the real field [12]. 

This implies that even in the absence of channel distortion and noise, and with perfect time and frequency 
synchronization, there will be some inter carrier (and/or intersymbol) interference at the output of the 
Analysis Filter Bank (AFB), which is purely imaginary (the notation in [9] is adopted): 

Y,s m ,S l )g P , q ^ = j{s) P m l' (2) 

and known as intrinsic interference [7]. Making the common assumption that the channel is (approximately) 
frequency flat at each subcarrier and constant over the duration of the prototype filter [9], which is true for 

practical values of L h and L g and for well time-localized g's, one can express the AFB output at the p th 

subcarrier and q th FBMC/OQAM symbol as: 



H d +jj Y H d (gY-'+r] 

p,q p,q J i m.« m»" \*» I m t n ' P>9 

m=0 n(m,ti)*(p,q) 



(3) 



Where H q is the Channel Frequency Response (CFR) at that FT point, and I and rj p are the 

associated interference and noise components, respectively, rj has been shown to be stationary (in fact, 

also Gaussian with zero mean and Variance (J 2 ) and correlated among adjacent subcarriers (see [5]). This 
correlation will be henceforth assumed negligible, a valid assumption for well frequency-localized prototype 
filters. 

A common assumption is that, with a well time frequency localized pulse, contributions to / only come 
from the first-order neighbourhood of (p,q) , namely = {(p±l,q±l),(p,q±l),(p±l,q)} . Moreover, the 
CFR is almost constant over this neighborhood, one can write (3) as 

Where , (5) 

p,q p,q J 111,11 \6 / m n p,q J p,q 

[»,«)eil M 

" P.I 

Hfi = HH(T,v)e J " i2lT " v - 2kv « + " ) dvdT (6) 
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Represents the channel coefficients at q' h symbol p th subcarrier frequency, and 

p-'i 

is the virtual transmitted symbol at (p , q), with 

u = y d (g) p -" ( 8 ) 

p,q Z-l ™,» \ 6 I m,n 

(m,n)efl i)ia 

being the imaginary part of the interference from the neigboring FT points. When known pilots are 

transmitted at that FT point and its neighbourhood ^ , the quantity in (5) can be approximated. This can 

be serve as a pseudo pilot [9] to compute an estimate of the CFR at the corresponding FT point, as, for 
example, 

"■« c "•" c 

p,q p,q 

This observation underlies the IAM schemes to be described in detail in the sequel, 
in. Iam Preambles With Guard Symbols 

From equation (9) that the preamble should be so organised as to result in pseudo-pilots of maximum 
magnitude [9]. Therefore, the training symbols surrounding the pilot d should be such that all terms in (8) 
have the same sign so that they add together. Moreover, this should occur for all frequencies p. 
To this end, one needs to know the interference weights (e) for the neighbors ( m, «)eQ of each 

X*- 7 / m,n ° v ' ' P*q 

FT point (p,q) of interest. These can be a priori computed based on the prototype filter g employed. 
Furthermore, it can be shown that, for any choice of g, these weights follow a specific pattern, which, for the 
definition of (p Q adopted here, and for all q, can be written as 

i-iys -p (-iys 

-(-D p r d pq (-iy r (io) 
i-iys p i-iys 

With the horizontal direction corresponding to time and the vertical one to frequency. The above quantities 
can be given by 

_ . 7* L t -1 L t -\ , 2ir 

P = e J M 2 £ g 2 {l)e' M (11) 



. M 



_ .2k V-l L g -l 
l=M 12 



(12) 



M , fe> (13) 
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and are positive and, clearly, smaller than one. Generally,/?,}/ >5. for example, y = 0.5 644 /? = 0.2393 
and S = 0.2058. for the g's employed in part V. This observation, along with the symmetry of the above 
pattern, will be seen to underly the development of the IAM method proposed in Section IV. 

A. IAM-R 

In order to simplify the task of generating pseudo-pilots of large magnitude, it was proposed in [15, 14] to 
place nulls at the first and third FBMC/OQAM symbols of the preamble, that is, 
a P o =a P 2 =0,p = 0,l,....,M-l. Then, the imaginary part of the pseudo-pilot at (p,l) will only come from the 
symbols at the positions ( p + 1, 1) • Obviously, these pilots should be OQAM symbols of maximum 
modulus, d. Moreover, in view of the pattern (10), they should satisfy the condition d +l t = - d l , for all 
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p, in order to yield pseudo pilots of maximum magnitude, namely \d , + j2Bd +11 | = d-Jl + 4B 2 . an 

example for M =8 and OQPSK modulation is shown in Fig 1(a). This method is called IAM-R to signify the 
fact that its pilot symbols are real. 

B. IAM-C 

Regardless, it can be verified from (5) that one can do better in maximizing the magnitude of the pseudo- 
pilots if imaginary pilot symbols jd , of maximum modulus d are also allowed in the preamble. The 
reason is that the corresponding pseudo-pilot would then be imaginary and hence of larger magnitude than in 
IAM-R, namely \d pl + u p J, provided that the signs of the pilot symbol and its neighbors are properly 

chosen so that d |W ( > • Such an IAM scheme was first suggested in [11], and was called IAM-I to 

signify the presence of imaginary pilots. The middle preamble vector comprises of triplets, with each of them 
following the above principle and the pilots in each triplet are otherwise choosen in random and 
independently of the other triplets. Hence, imaginary pseudo-pilots result in only one third of the subcarriers, 
whereas the rest of them deliver complex pseudo-pilots of smaller magnitude, d 1(1 + B)+ j/3\. With a 

minor but significant modification, one can improve upon this and obtain pseudo-pilots that are either purely 
real or imaginary at all the subcarriers. The idea is to simply set the middle FBMC/OQAM symbol equal to 
that in IAM-R but with the pilots at the odd subcarriers multiplied by j, as shown in [8] and independently in 
[6]. All of the resulting pseudo-pilots are then of maximum magnitude, d |l + 2/9 1. Again, for the M = 8 

OQPSK example, this will result in the preamble shown in Fig. 1(b). This scheme was called IAM-C. Note 
that this preamble is, firmly speaking, not OQAM. Regardless, it can still be fed to the synthesis filter bank 
and perfectly reconstructed at the analysis bank. 

C. IAM-C EXTENDED(E IAM-C) 

It was shown in [6] that the pseudo-pilots generated by IAM-C are of the maximum possible magnitude. 
However, this optimality is constrained to the 3 -symbol preambles with nulls at their sides. The property 
mentioned at the end of part II suggests that an alternative 
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Fig 1 : preamble structures for (a) IAM-R (b) IAM-C and (c) E IAM-C 

Preamble structure which employs the side symbols as well could yield pseudo pilots that are stronger than 
those of IAM-C. It turns out that such a preamble can really be formed in an analogous way as previously, by 
using the left- and right-hand sides of the neighborhood (10) as well. Thus, at an odd-indexed subcarrier p, 
with the pilot dj in the middle, one should place +d at the right hand side, and its negative, +d at the left 
hand side. These will then contribute idyj each to the interference. In an analogous way, at an even-indexed 
subcarrier p, the middle pilot, say +d , will get a total interference of +dj from its (p,0) and (p,2) 
neighbours, if these are chosen as + dj and + dj ■, respectively. It can then be verified that, because of (10), 
the interference components from the corner neighbors cancel each other, and hence have no contribution to 
the pseudo-pilot. However, this loss is not significant since d is significantly smaller than g. The resulting 
pseudo-pilots are again either purely real or imaginary, with magnitude c/|l + 2(/3 + /)|, which is clearly larger 
than that achieved by IAM-C. 
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IV. Optimal Iam-C(0 Iam-C) 

Inspired by IAM-C and observing that with the triplet using «(1 + <J)± jS nave smaller power 
compared with djj , one may think of finding the optimal complex valued d m , , carefully so that the 
demodulated symbols di"] will have the maximum power. 

To further increase the power of demodulated preamble symbols djf} , the two neighbouring columns can 

further more be utilized by assigning proper complex valued symbols rather than zeros. For convenience, the 
following notations are used to simplify our derivation. 

A g (0,±v ) = S, A g (±T ,0) = p, A g (±t ,±v ) = r , A g (±T ,±2v ) = n 

For distinct pulse shapes g(t) the above parameters are distinct in general. For well designed pulse shapes, 
A (nr ,mv ) attenuates very fast with increasing m, n. For the demodulated symbol a£j , a weighting 
matrix W can be formulated to calculated prominent interference elements from neighbouring symbols. As 
shown in (7), the sign of the elements in wwill be distinct for distinct value of / . Therefore two different 
realisations are listed as W, „ when / is an even number and W... . 
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w kven an< ^ w hdd on ^y differs from each other by altering the sign of all the elements in the second row and 
the fourth row. Since the first and the third column of the weighting matrix w also furthermore contribute to 

■tic) 

intrinsic interference, more power will be conveyed to d k , by assigning proper complex valued symbols 

rather than zeros. The symmetric and fast decay nature of the weighting matrix w in the central column are 
repetition of a quadruplet, the symbols in the other two columns should also be grouped into quadruplet. As 
shown in Appendix A, given the weighting matrix w , the optimal preamble sequence is as shown in Fig 2 
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Fig 2: optimal preamble structures for (a) optimal-IAM \J even ) (b) optimal-IAM (/ , . ) 

This new method was named as optimal-IAM and the preamble structures are shown in Fig 2, and the 
preamble structure appears repeatably to formulate the preamble sequence. 

V. Simulation Results 

where N =512, 1024, 2048 FFT/IFFT was used with a convolutional channel coding at rate 1/2, we simply 
use a relatively small FFT/IFFT .The first three IAM preambles, namely IAM-R, IAM-I and E- IAM-C, have 
the same energy on preamble sequences and share a similar structure (0 X 0) and there will be frequently 
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compared with each other in our following discussion. The last IAM preambles, namely Optimal- IAM-C, 
have the same structure as defined by O-IAM-C (J even ) and O-IAM-C (/ orfrf )but three times the energy of the 
other three preambles. The normalized MSE \\ H _ H - || ! where His the channel frequency response vector & 

II" If 

H its estimate, is plotted with respect to the E b I N a (SNR). BER simulation versus e, I N (SNR), with 

E, the useful bit energy and N Q the mono-lateral noise density, has been carried out for different preamble 

sequence, as shown in Fig.l&2. The O-IAM-C method can always outperform E-IAM-C by at least 1.24 dB. 
Fig. 4 shows the MSE performance of the methods under study for M= 512, 1024 & 2048 subcarriers with 
an overlapping factor of K = 4 are considered in Figs. 4(a), (b) and (c), respectively. Observe that in three 
cases O-IAM-C outperforms the other schemes in the whole SNR range considered. Fig. 3 shows the BER 
performance of the methods under study for M=512, 1024 & 2048 subcarriers with an overlapping factor of 
K=4 are considered in Figs. 3(a), (b) and (c), respectively. Observed that in all the cases O-IAM-C 
outperforms E-IAM-C by at least 1 .24 dB gain. 



VI. Conclusion 

The problem of preamble -based CE in FBMC/OQAM systems was revisited in this paper, with a aim on the 
interference approximation idea. The three most well renowned and efficient IAM variants, IAM-R, IAM-C 
and E-IAM-C, were reconsidered along with their associated preamble structures. Emphasis was put on IAM- 
C and its optimality within the class of 3- symbol preambles with all zeros side FBMC symbols. The possible 
gains from resting the last mentioned constraints of all zeros at the sides were investigated, on the basis of 
fully exploiting the symmetries in the interference weights of the first-order FT neighbours. An innovative 
preamble structure of this kind, optimal -IAM-C, was developed as an extension of the E- IAM-C idea and 
shown to result in pseudo-pilots of significantly larger magnitude. The superiority of O-IAM-C was 
demonstrated via results of simulation in a realistic (including interference from data) and fair scenario. 

Appendix A 

A. Derivation of the optimal IAM preamble structure 
For odd 1, given the preamble matrix 



x, 1 

*2 j 



1 X - 



J X > 



(15) 



Where X ; < l,i = 1, ,8 are complex valued variables, the demodulated symbols at corresponding points 

are as follows: 

a\ c) = 1 + 28 + j [f3(x 5 - x 1 ) + y(x 2 + x i + x 6 + x g ) + 2r](x 7 - x 3 )] 

>'i 

af = + 28 + P(x 6 - x 2 )+ r(x l + x, + x s + x 7 ) + 2r}{x t - x 4 )] 



= -1 - 28 + j [/3(x 7 - x 3 ) + y{x 2 + x 4 + x 6 + x g ) + 2rj{x s - Xj)] 



(16) 



-j = ~ 28 + p (x g -x 4 ) + y(x, + x 3 + x 5 +x 1 ) + 2ri(x fi -x 2 )] 
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Then the preamble design problem becomes a optimisation process which try to find a realisation of 
x. , i = 1 8 that can maximise the minimum power of the demodulated symbols, i.e. , that can maximise 

the minimum power of the demodulated symbols, i.e, 

{x i ,/ = l,...,8} = arg max g min{|a 1 ( °| 2 ,|aj°| 2 ,|a k 1 ) | 2 ,|a k y ) | 2 } (17) 
Under the constraint that |x ( . I < 1, i = 1, , 8 . 

In order to maximise the power, all the elements in each a ' in (16) should be coherently added, that is, 

phase(x 5 - x,) = phase(x 2 + x 4 + x 6 + x 8 ) = phase(x 7 - x,) 

phase(x 6 — x 2 ) = phase(x { + x 3 + x 5 + x 7 ) = phase(x s — x 4 ) (18) 

Given (16), it is easy to figure out that the condition that maximise | 2 (e.g. when jy l >0) will 

minimise | fl <^ | 2 , and vice versa. The same situation happens for | a < c) | 2 and \a {c) . | 2 . Therefore the optimal 

solution to maximise the minimum power is to set y l ,y i purely real valued and y 2 , y 4 purely imaginary, 
that is, 

phase(y l ) = phase(y 3 ) = on 

phase(y 2 ) = phase(y 4 ) = n I 2 or -nil 

Therefore following equations can be formulated as 

{3(x 5 ) = 3(x,), 3(x 7 ) = 3(x 3 ) y t ,y 3 purely real 
*R(x s )5?(x,) < 0,5?(x 7 )9?(x 3 ) < coherent addition 

*R(x,+x 3 +x 5 +x 7 ) = ( 19 ) 
9J(x 6 ) = fft(x 2 ),fft(x 8 ) = 9?(x 4 ) y 2 ,y 4 purely imaginary 
0(x 6 )3(x 2 ) < 0, 3(x 8 )3(x 4 ) < coherent addition 
3(x 2 + x 4 + x 6 + x g ) = 

Where S J?(.) denotes the real part operator and 3(.) denotes the imaginary part operator. It is clear from 

(17), (18) and (19) that the power of y.J = 1,....,4 can reach the maximum only if |x ; | = l,i = 1, ,8 , 

i.e., each element in the preamble utilises the maximum power available. Hence (5.21) can be rewritten as 

3(x 5 ) = 3(x,), 5R(x 5 ) = -9?(xj) => x, = -x* 

3(x 7 ) = 3(x 3 ) 9?(x 7 ) = -9?(x 3 ) ^ x, = -x* (20) 

< 

9?(x 6 ) = «(x 2 ), 3(x 6 ) = -3(x 2 ) ^ x 2 = x* 

9?(x g ) = 9?(x 4 ), 3(x g ) = -3(x 4 ) x 4 = x* 

Hence the preamble matrix A lodd described in (15), we get is 
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Fig3: BER vs. SNR with QPSK modulation .Rayleigh channels were assumed, with 
(a) Filter banks with M=512,K=4(overlapping factor), (b) Filter banks with M=1024, K=4 (overlapping factor) and (c) Filter banks with 

M=2048, K=4(overlapping factor) 
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Fig4: Estimation (NMSE) performance of preamble based methods .Rayleigh channels were assumed, with 
Filter banks with M=512, K=4( overlapping factor), (b) Filter banks with M=1024, K=4( overlapping factor) and (c) Filter 
banks with M=2048, KM^overlapping factor) 
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